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ABSTRACT 

Using many-body techniques we obtain the time- dependent Gaussian approximation for 
interacting fermion-scalar field models. This method is applied to an uniform system of 
relativistic spin-1/2 fermion field coupled, through a Yukawa term, to a scalar field in 3+1 
dimensions, the so-called quantum scalar plasma model. The renormalization for the resulting 
Gaussian mean-field equations, both static and dynamical, are examined and initial conditions 
discussed. We also investigate solutions for the gap equation and show that the energy density 
has a single minimum. 
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I. Introduction 



The self-consistent mean-field method remains being one of the few analytical tools avail- 
able to investigate a variety of problems in quantum field theory which cannot be handled in 
perturbative theory. These include static problems in the context of spontaneous symmetry 
breaking as well as time-dependent problems involving the dynamics of the inflation- 
driving scalar field in the early Universe and the dynamics of relaxation process in the 
ultra-relativistic heavy-ion experiments 0. In its most used form the approximation is imple- 
mented, in the case of boson fields, through of a time-dependent or independent variational 
principle using a Gaussian trial wavefunction in the Schroedinger picture 0]. However, sys- 
tematic corrections to this approximation are still an open problem. Furthermore, for the 
case of fermion fields, the trial wavefunction is not straightforward to take in a Gaussian form 
i- 

On the other hand, the mean-field description of systems of interest can be obtained from 
the Heisenberg equation 

t{d)=Tr\0,H]ro , (1) 

where (9 is a one-body operator and H is the Hamiltonian of the system. This can be kept 
under direct control when one takes a gaussian-like density. In a formulation appropriate for 
the many-body theory we have 



exp 






Tr |exp 




} 



•^0 — 7 f 1 r~f TT ■ {^) 

Tr |exp AjViVj + BijvlVj + j 

In d) the 7]} {7]i) are boson or fermion creation (annhihilation) operators of a particle in 
orbital i and J-'o can be, of course, rewritten in a diagonal form when one uses the Bogoliubov 
quasiparticle operators (cf. Eqs.(^ and (|l3)). It is not difficult to see that this scheme 
corresponds to the usual gaussian variational approximation in the functional Schroedinger 
picture. There are, however, some important differences. For instance, the usual technique of 
nonrelativistic many- fermion problem can be readily extended to the Dirac field theory 0. 



Thus, instead of using anticommuting Grassmann variables to write Gaussian wavefunctional 
we describe the pairing correlations by a suitable Bogoliubov unitary transformation. 
Furthermore, this could be the starting point for the further approximation scheme beyond 



the Gaussian [13|. In recent publications we have applied this scheme to the 0^ theory 
and Chiral Gross-Neveu Model (hereafter referred to as I, II). 

The main point of this paper is to apply this method within the context of interacting 
fermion-boson models. As a first step towards this end, Takano Natti and de Toledo Piza 
0] have obtained relevant dynamics for the Jaynes-Cummings Hamiltonian, a well known 
model in quantum optics [|l^. This can be seen as + 1 dimensional quantum field theory 
known as relativistic scalar plasma 0] . Exact numerical results have been useful in this case 



to assess the quality of various approximation schemes. In particular, we have shown that 
the method described here is a good approximation whenever a moving Gaussian has enough 
freedom to mimic the complicated evolution of the exact wave function. In this paper, we 
report an aplication of the same technique to describe the real-time evolution for a fermion 
field interacting with scalar boson in 3 + 1 dimensions. An outline of the paper is as follows. 
In Sec. II we extend the discussions of I and II to a system countaing interacting boson 
and fermion fields. Sec. Ill will illustre this scheme in the simplest context of scalar plasma 
model. In Sec. IV we discuss the self-consistent renormalization for the resulting equations 
in the equilibrium situation. The renormalization of the time-dependent equations is then 
considered in Sec V and finite initial conditions discussed. 

II. Kinetic of Interacting Fermion-Scalar Systems 

The present discussion for the real-time evolution of theories containing both scalar and 
fermion fields basically follows the earlier works in the context of 0^ and Chiral- Gross-Neveu 
models. The basic idea of our approach is to focus on the time evolution of the expectation 
values of linear, and bilinear field operators such as , il){x)ip{x) , ip{x)ip{x) and 

etc, referred to as Gaussian observables. Many technical details of the method can be found 



in I and II, and we restrict ourselves here to a few key steps of the derivation of equations of 
motion for the present case. 

Let us consider first the bosonic sector of the system. The Heisenberg field operators (f){x) 
and n(x) are expanded as 

1 



(3) 



n(x,t) 



E 



where &p(t) and bp{t) are the usual boson creation and annihilation operators satisfying the 
standard commutation relations at equal times 

h{t),bl,{t')U> = Vp' , [bl{t),bl,it')U, = [6p(t), V(^')]*=*' = . (4) 
In Eq.(|^) V is the volume of the periodicity box, 

(Po)^ = + fi^ and px = pot — p.x , 

where Q is the mass parameter and will be fixed later. 

In order to deal with condensate and pairing dynamics of the scalar field we follow Sec. 
II of I and define a unitary Bogoliubov transformation as 



cosh Kp — sinh Kp — 



sinh Kp + cosh Kp + 





dp 




. dip _ 



(5) 



where dp is the shifted boson operator 



dp{t) = bp{t) - B{t)6p^o with Bp{t) = {bp{t)) = Tr^^ [bp{t)T] 



(6) 



Here and in what follows the symbol Trgp denotes a trace over both bosonic and fermionic 
variables. Partial traces over bosonic or fermionic variables will be written as Tr^ and Trp 
respectively. The Bogoliubov coefficients are determined from the secular problem given in 



(2.8)-(2.10) of I to yield {PpP-p) = 0. The resulting eigenvalues, z/p = {P^Pp), are quasi- 
particle occupation numbers. 

For the fermion field, we follow the discussion of II and decompose the field operators in 
Fourier components. 



M\^/2 1 



(7) 



(2) /+\„ik.x 



where Oks^*^) ^^'^ '^ks(^) (^kii^)] fermion creation and annihilation oper- 

ators associated with positive [negative] -energy solutions ui(k, s) [u2(k, s)] and satisfy the 
anticommutation rules 

{4f W,«L''kO}t=t' = S^M'^s,s'Sx,x' for A, A' =1,2 



In Eq.(|^) the spinors are normalized to ko/M and 



{koY = k^ + and kx = kot - k.x 



(9) 



where M is a mass parameter for the fermions. 

Adapting Eqs. (3)-(ll) and (30) of II to the present case of 3+1 dimensional field we can 
write the Bogoliubov quasi-particle operators as 



a 



(2) 
k,s 



a 



(1) t 



^k,s 



a 



(2) T 







COS(/?k 




-e^'^^ sin 







cos(^k — e sin 



g^k sin yjj^ cos (/9k 



e sin 






cos v^k 



.(1) 



,(2) 



,(1) 



,(2) 



(10) 
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and require (a^kjQJk's) = in order to determine this unitary transformation. The quasi- 
fermion occupation number will be given by u^^l = (a^^]^a^^]). In ([T0|) we used \s) = 
(— )^/^^^| — s) to write the spin labels. 

The next step is to obtain the equations of motion for the condensate, (0) and (11), and for 
the Bogoliubov parameters given in (|D-(|D and (|10|). This proceeds by using the Heisenberg 
equation of motion (|l|), where O are performed by linear or bilinear combinations of creation 
and annihilation parts of the field operators [see (2.17) of I and (5) of II for details]. For the 
bosonic variables we have 

= Tr^J6.H]Tn (11) 



2(11) = Tr,,[U,H]J^o (12) 



z>p = Tr,Mp,H]J^o (13) 



Analogously, the dynamical equations for fermionic variables read as 



Tr.A[<l<l^H]:Fo) (15) 



Tr BF ([a'll s<^kl^ H]J^o) 
[«0k + 7k sin cos v^k] e = — '—t^. . (16) 

1 - ^k,s - 

Our implementation of the Gaussian mean-field approximation consists in assuming a 
factorized form for the density matrix Q 

^o(t) = ^o'-^o' • (17) 
In Eq. ([l7|) the subsystem densities JFq and JFq are gaussian densities given in Eq.(H), now 
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rewritten in the diagonal forms using the Bogohubov quasiboson and quasifermion operators, 



^o=U-r^iT^] . (18) 



k,s,A 



Substituting into (|T^) and (|1^) one verifies immediately that the occupation numbers are 
constant, 

z>p = rrBp[/?J/?p, H]J^o = ^>k,A = TrBp[4!]^«M, = . (20) 

Thus, Eqs.([TT| ) -([16| ) form a set of self-consistent equations of motion to the Gaussian variables. 
Its implementation for a specific model consists essentially in taking the traces of appropriate 
commutators in the Fock space. The fact that occupancies are constant implies that the 
entropy function associated with JFq, i.e., 5* = — Tr^^pJ^o log J-'o, does not change in time. One 
recovers therefore the usual isoentropic character of the Gaussian approximation. 

A method to improve the mean-field approximation in field theory was discussed in I. The 
approach follows the line of thinking of a time-dependent projection technique proposed some 
time ago by Willis and Picard in the context of master equation for coupled systems and 



extended later by Nemes and de Toledo Piza to study nonrelativistic many-fermion dynamics 
ITSl . The method consists essentially in writting the correlation information of the full density 
of the system in terms of a memory kernel acting on the uncorrelated density Fq, with the help 
of a time- dependent projector (see Sec. Ill and Appendix A of I for details). At this point, 
an systematic mean-field expansion for two-point correlations can be perfomed. The lowest 
order reported in this work corresponds to the results of the usual Gaussian approximation. 
The higher orders describe the dynamical correlation effects between the subsystems and are 
expressed through suitable memorial integrals added to the mean-field dynamical equations. 
Thus, the resulting equations acquire the structure of kinetic equations, with the memory 
integrals performing as collisional dynamics terms which eliminate the isoentropic constraint. 



This higher order approximation has been implemented for the case of + 1 dimensional field 
model and will not be done in this paper. 

III. Effective Dynamics of Relativistic Scalar Plasma 

The procedure discussed in the preceding section is quite general and can be, in prin- 
ciple, utilized for systems containing scalar and spin-1/2 particles. This is not, however a 
straightforward task, since divergences appear in the equations require a discussion of the 
renormalizability of specific models of interest. In Ref. jl^ properties and stability of possi- 
ble vacuum states for several fermion-scalar models utilizing the Gaussian-effective-potential 
(GEP) approach have been studied. It is shown that in 3+1 dimensions the models seem 
doomed to instability and the GEP is bounded below only if the Yukawa coupling vanishes. 
We have reached a similar conclusion in an early stage of this work. 

In order to give a well defined physical application, we shall adopt here the relativistic 
scalar plasma model widely used within the context of ultrahigh dense matter [|I^]. This 
system consists of relativistic fermion gas interacting through a massive scalar meson only 
and corresponds to one of the simplest quantum-field theoretical models used to discuss 
the relativistic dense matter in the contexts of heavy-ion collisions and the high-density 



astrophysical system |]1T|. In this application, one considers the mean-field contributions, 
(0), are much larger than the fiutuactions, (0^) — (0)^- Therefore, the bosonic pairing is 
ignored and Bogoliubov transformation reduces to identity. We shall show in the next 
sections that the resulting equations of interest in this case can be made finite by introducing 
appropriate counterterms. 

We now proceed to apply the technique discussed in Sec. II to the scalar plasma model. 
The dynamics is governed by the Hamiltonian (we use the notation: = / d^x) 



H 



I n, 

(21) 



n = -^(2^.^-m)7/'-c/7/'0v^ + — f(47rn)2 + |a0P + /iVl 



where the parameters m and /i are, respectively, the mass of fermion and scalar particles 
and g is the coupling constant. The last term of this expression contains the counterterms 
necessary to remove the infinities occurring later ||15|| , 



A ^ 5fx\2 C . D 



(22) 



where the coeffcients A, Sfi'^, C and D are infinite constants to be defined later. 



Introducing these ingredients into Eqs.(|ri|), (|l^) and ([iq ) one easily obtains dynamical 
equation for the condensate 



47r(n) 



(23) 



(ri) 



1 

An 



A + ^GiSl) 



1 



I" 1 

/ — [M COS 2v9k + k sin 2v9k cos 7k] (1 - t'?] - 4' 



(1) „ „(2)x 

s J 1 



(24) 



where k = |k|, and two real equations for the fermion pairing 

k 



ko 



(M — m) sin 7k 



(25) 



2(k2 + mM) 2(M-m), 
sm2(y9k7k = ; sm2(^kH ] A; cos 2(^k cos 7k 



In 



ko 



we have introduced the notation 



l + 2z/„ 



p 2^/p^TW 87r2 I P 



ko 



y/eQ J Jp po 



(26) 



(27) 



where Ap is a momentum cutoff used in evaluating the integral. In (i^)-(i3) m = m — g{(f)) 
stands for the effective fermion mass. 
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Another physical quantity of interest is the energy density of the system, 
{H) 1 



V 



V 



-TrHTn 



E 



(k^ + rfiM) im- M)^ 

— cos 2y9k H ; k sin 2v9k cos 7k 



ko 



'1 (1) (2)n 



+- 



in 



47r 



(0) 



Stt 



(5/1^ + —Gin) 



247r 



D 



1 



9o7r Stt 



f28) 



An important feature of this scheme is that the mean energy is conserved, a property which 
can be verified exphcitly using the equations of motion p3|)-(p6D. Notice also that the results 
above contain divergent integrals, so that the coefficients A, 6fi'^, C and D have to be adjusted 
appropriately in order to give a well defined dynamics. We shall focus on this point in the 
next sections. 



IV— Static Equations and Renormalization 



This section considers Eqs. (PB|)-(P^) in the equilibrium situation. We shall investigate 
the solution of these equations and study renormalization conditions. Hence, we set 



7k = 0k = (0) = (n) = 



(29) 



Since our main goal here is to set up a Gaussian mean field approximation for interacting 
fermion-boson fields, we shall, for simplicity, not consider the matter contribution in the 
present calculation. This can be done in the case of Gaussian mean-field approximation 
by assuming an appropriate statistical distribution to the occupation numbers z/p and v^l 
as initial conditions in (|l^) and (|T5l). Thus, for the rest of discussion we have vanishing 
occupancies, z/p = u^l = 0, and the sum in spin gives a numerical factor, J2s = 2. With these 
assumptions, Eqs. (|23|) -(p6D become 
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(n)L, = o, 



(30) 



+ 



+87r5f / -3- (Mcos2v9k|cq + /csin2v9k|cqCOS7k|eq) = , 



(31) 



(M - m) sin7k|eq = 



(32) 



(k^ + mM) sin2v9k|cq + (M - m) cos2v9k|cq cos7k|eq = 



(33) 



where fh is now the effective fermion mass calculated at (0)|cq- Solutions of this set of 
equations correspond to a stationary points of the energy density, and include the ground 
state of the system in this approximation. 

A. Equilibrium conditions and mass parameters 

Considering equations (^2|) and ( PB| ) one sees that two situations have to be analysed. 




Using this solution in Eq.(P5D we have 

(k^ + rhM) 

cos2v?k|cq = ± 



sin2y9k| 



M -rfi 



(34) 



where ko = Vk^ + m^. Substituting this result into ( pi]) one finds for this case the following 
gap equation 



+ 



(0)|eq + ^{<P)t + §i<l>)t = ±lQ^9mG{m) , (35) 



where G{m) is the divergent integral defined in (]27|) with = 0. 
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ii) M = m 

In this case 7k can have any value and (^3D has as solution 



sin2v9k|eq = , cos2v9k|cq = ±1 (36) 



Comparing now (Q) with (|36[) one notices immediately that the former includes the second 
solution as a particular case, where the mass parameter M is defined here by the equilibrium 
condition and Bogoliubov transformation (|1^) reduces to an identity matrix. 

The above discussion shows that the equilibrium value (0)|cq of the condensate, or the 
effective fermion mass m are independent of the initial choice of the mass paramete M. In 
fact, the effects of pairing transformation will be such that the effective fermion mass is always 
given by the equilibrium value calculated with (|35|). (See also (^) below for the mean-field 



energy). Henceforth, we shall thus take M = m without loss of generality. 

B. Counterterms 

The arbitrary parameters A, 6fi, C and D can be fixed easily by adjusting the coefficients 
of (0)|oq in both of sides of (|35|) to give a finite gap equation. Thus, taking account the ± 
signs in ( p5D we choose the following self-consistent renormalization prescription 

D = ±487r^U(m) , (37) 

6^^ = T24n^g^L{m)G{iJ,) T lQng'^G{0) ± 24nm^g'^L{m) , (38) 

C = T^8nmg^L{m) , (39) 

A = ±24:7img^L{m)G{fi) ± 167imgG{m) , (40) 

where 



, . /• 1 1 2Ap ^ ^ 

^("^^ = i 2k^(k^ + m^)V2 = 4^^ ^ • (41) 

Furthermore, fi and m are the mass scales for boson and fermion fields respectively. 
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Substitution of these counterterms into (^) produce the appropriate cancelations which 
render the equation finite. A combination of type L(m)[G(/i) — G'(fi)] comes from the first 
two terms. Since fl is an arbitrary expansion mass parameter, one can remove this divergence 
by setting Q = fi. The resuhing finite gap equation is 



m/ 2 



. (42) 



The above equation together with (p^)-P^) determine the stationary points of the model in 
the Gaussian mean-field approximation. 

C. Mean energy and stationary solutions 

We examine next the energy density when it is stationary with respect to the fermion 
variables. With the help of (p4|), Eq. f^Sf ) becomes 



(El 

V 



(<^k|cq,7k|cq, (0)) = T 



{4>) + 



in 



(43) 



247r'"' 967r 

Substituting next the counterterms (p7D-(^0D into this result, we have the renormalized version 
of mean-field energy as function of mean-field value (0) or fh 



V 87r2 



vr/i 



(m — mY + In 



+ 



m 



(44) 



where we have added appropriate constant in order to get (iJ)((0) = 0) = 0. One can now 
discuss the possible solution of (|^) by analysing the minima of (^31) . 

Let us define x = g{(f))\,Jm and E{x) = {Sn"^ /m^){H) /V. Eqs. (H) and (||) can then 
be written respectively as 
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(45) 



and 




^ ln(l-x) + 



4 



n 



4 



1 



(46) 



The combination "^^^ has been used as an effective couphng constant in Ref . jlSl . 

The behavior of E{x) is shown in Figs. (1)- (4) for several combination of fi/m and g^. 
Notice first that this function has the domain at < x < 1, which is the physical range 
for the fermion mass. The point x = corresponds to (0) = or m = m and a; = 1 is 
the case when the effective fermion mass m = 0. Qualitatively, the results indicate that the 
system always presents a single minimum. Figs.(l)-(2) show E{x) for several values of 
with /i/m fixed (see figure captions for numerical values of the parameters). The positions of 
Xmin indicate that the vacuum approaches x = when we decrease g"^. In the limit of g^ — > 0, 
one gets Xmin —>■ 0, as shown in Fig. (3). In this case, m ^ rh is the optimal fermion mass, 
as it must be in the free field theory. On the other hand, Figs.(3)-(4) plot the function E{x) 
keeping the value of the Yukawa coupling g^ fixed, but with different values of the ratio fi/m. 
Comparing these two curves one sees Xmm ^ when fi/m —>■ oo. In other words, when the 
meson mass is large, the force range is small, as usual in the Yukawa theory. In the limit of 
infinity /i, the fermion particles of the system cannot interact. It can be seen also from the 
AH's formula, where m//i plays roles of effective coupling constant. The above discussions 
suggest that the field has always a stable vacuum. This means that there is a finite range 
around the minimum where the dynamics of the system is well defined. In the next section, 
we shall therefore discuss renormalization conditions for the time-dependent equations. 

V. Renormalization and Initial Conditions for the Time-Dependent Equations 

The last section has discussed the problem of renormalization for the vacuum sector of the 
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relativistic scalar plasma model in the Gaussian mean-field approximation. We have shown 
that the physical quantities can be made finite with the counterterms introduced in (p2D and 
Here we shall consider an off-equilibrium situation and study renormalizability for 



the Gaussian equations of motion for this model. (See, e.g., Refs. |16, 17] for the issue of 
renormalization of time-dependent equations in 0^ field theory) 

A. Equations of Motion and Initial State 

In order to pose the problem more clearly, let us begin by rewriting Eq. (^) , with the help 
of (^-(0), as 

(0) = -/i^(0)-8vr^^(t) . (47) 

where 

nt) = [ -^k(t) = / 

7k Jk 



m „ k . ^ mit) m^it) 

— cos 2y9k + -r~ ^^'^ -^V^k cos 7k + 



(48) 



ko ko k 2k'^ko 

Notice that the integral is divergent unless JF^ decreases faster than k~^ . For simplicity, we 
shall ignore here the possibilities of fractional powers, logarithmic or oscillatory behavior and 
take the convergence condition for the integral to be 

limJk = cfA) ■ (49) 



Thus, the allowed domains for the dynamical variables in momentum space suffer strong 
constraints, and our investigation of renormalizability of the equations of motion consists in 
analysing the large momentum behavior of JF^ from this point of view. In order to satisfy 
(^) at all times the pairing variables and 7k must have their time evolution restricted. 
These restrictions are analysed better by making the variable changes 



/cq cos 2y9k = —mR]^ — k\J kl — (50) 

or, equivalent ly. 



kl sin 2v9k = -/c-Rk + m^J kl - Rl , (51) 
15 



and 



cos 7k = 1 — VFk 



(52) 



In terms of these new variables we have 



k ' 



(53) 



where 



7r(l) 



i?k rh{t) rfr-{t) 
ko k 2k'^ko 



(54) 



(55) 



We see in J^^^ that the leading term of -Rk^oo must be m in order to cancel the quadratic 



divergence. Thus, we define 



lim R]^ = m + Sy^ 

fc— >oo 



(56) 



valid for large k. Using this in (^) one gets 

,2 



lim J"^^^ - - — 

fe^oo ^ kk^ik + ko) 2k'^kQ k^ 



/co L k{k + fco) 

There is still a logarithmic divergence left. To keep this under control we must have 



(57) 



lim S'k , , , , , 



(58) 



in order to get finite /k^k^'*- 



We examine next the convergence condition for /k^k • Using (pO) and (BSf) in (B3) we 



find the asymptotic expression for J^^^ 



lim JF, 



(2) 



k^{w? — rn^) Q ( ^ 
kQ(km + kom) \k^ 



(59) 
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This will constrain the large-/c behavior of W]^ to be 



(60) 



in order to satisfy (|^). 

The asymptotic conditions (^), ( ^8]) and ( |60| ) restrict the dynamics of the Gaussian 
variables in the moment space. Whenever the Fourier spectrum JF]^ has such properties 
Eqs.(|47|)-(^8D are well defined, and one ought to use these also as a criterion for choosing initial 
conditions. Since, however, we are dealing with a time dependent problem, the question we 
want to address to ourselves is whether the nonlinear evolution will distort the asymptotic 
behavior built into the initial state. From the energy viewpoint, we argue that, if the relations 
(p6|), and (pOf) result a finite energy density, its conservation will enforce these relations 



at all time. In fact, as we shall show below, the asymptotic behavior obtained in (|56D, (|58| ) 
and (IBO) are suffcient to make fBB) finite. 



B. Asymptotic Conditions and Energy Density 

In terms of new variables i?k and PFk the energy density can be rewritten as 



V 



where 



(1) 



1 



61T 



2k 



Jk Jk 



(2) 



Rl + 



ko 



— 2 

m 



- A 

m 



k Ak^ko 



(61) 



(62) 



(2) 



{m — m)k 



1.3 



W^{-kR^ + mJki-Rl 



Substracting an unimportant constant 2k from we find, after a little algebra 



4" = -k 



ko k 



2 k rm^ 
ko 



ko{ko + Jkl - RIY 



from which it follows that 



lim £] 



(1) 



mm 



kl \ k{k + ko)^^^ ^ ko 



m 



m / 1 
h C — 

4A;3 Akl W 



(63) 



(64) 



(65) 
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where we have used (|56D . It is clear now that this expression decrease as fast as k ^ and 
/]^ S^^ is finite. Comparing next (^3]) with (|55|) one notes immediately that 

lim ~ lim J^l^l (66) 

fc— >oo fc— >oo 

Hence, Sj^^ is finite when the condition (|^) is satisfied. 

The above qualitative discussion shows that the prescriptions given in (|56D , ( ^Sf ) and 



(|60|) can make energy density finite. They will thus define the initial states and the system 
will develop a well defined time evolution within a limited moment space. Because of the 
energy conservation, the asymptotic behavior of the Gaussian variables essentially will not 
be modified by the dynamics, otherwise this will cost infinite amount of energy from the 
system. This discussion has not considered, however, possibilities of eventual runaway, caused 
by eventual unbounded below behavior of the potential. The studies of this problem will 
require further detailed analysis of (H) as multidimensional function of m, i?k and W]^, (one 
coordinate for each k), and verification of its funcional properties [|T^. 

VI. Concluding Remarks 

In summary, we have presented in this paper a framework to treat the initial- value problem 
for interacting fermion- scalar field models. The method allows one to describe the real-time 
evolution of the fields in terms of the dynamics of few observables yielding a set of self- 
consistent equations for expectation values of linear and bilinear field operators. Although 
the procedure is quite general, we have however implemented the calculation within the 
simplest context of relativistic scalar plasma system. We have shown in detail that the usual 
form of renormalization also applies to the present nonperturbative calculation and we have 
obtained finite expression for energy density. A simple numerical calculation suggests that the 
system has always a single stable minimum, although further investigation will be necessary 
for other oscillation modes. The standard approach to this question is throught the use of 
RPA analysis, where the stability is indicated by its the eigenvalues. It is interesting to 
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mention here that the excitation modes described by the RPA equations are the quantum 
particles of the field. In fact, the physics of one meson and two spin-1/2 fermion can be 
investigated from this equation. We have also discussed the renormalization for the time- 
dependent equations. Using energy conservation as the key, we found that there are a finite 
range around the vacuum where the dynamics of the system is well defined. 

Finally, we comment that systematic corrections to include dynamical correlation effects 
to the present mean-field calculations can in principle be readily applied with the help of 
a projection technique discussed in I and Ref. 0. In this case, the occupation numbers are 
no longer constant and will affect the effective dynamics of the Gaussian observables. The 
framework presented here serves also as groundwork to finite density and finite temperature 
discussions [|^. In particular, a finite- matter density calculation beyond the mean-field 



approximation allows one to study collisional observables such as transport coefficients ||19|| . 
The extension of this procedure to explore nonuniform systems is straighforward but lengthy. 
In this case, the spatial dependence of the field are expanded in natural orbitals of extented 
one-body density. These orbitals can be given in terms of a momentum expansion through 



the use of a more general Bogoliubov transformation |20 
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Figure Captions 

FIGURE 1. The behavior of the ground-state mean-field energy density E{x) of the scalar 
plasma system as a function of fermionic effective mass x — g{(f)) /m — 1 —m/m for any values 
of the coupling constant g and mass scale fi/m = O.l fixed. 

FIGURE 2. The behavior of the ground-state mean-field energy density E{x) of the scalar 
plasma system as a function of fermionic effective mass x — g{(f)) /m — 1 — m/m for any values 
of the coupling constant g and mass scale ix/m — 2 fixed. 

FIGURE 3. The behavior of the ground-state mean-field energy density E{x) of the scalar 
plasma system as a function of fermionic effective mass x — g{4>) /m — 1 — m/m for coupling 
constant g"^ — tt/IOO and mass scale ii/m — 2. 

FIGURE 4. The behavior of the ground-state mean-field energy density E{x) of the scalar 
plasma system as a function of fermionic effective mass x = g{4>)/m = 1 — m/m for coupling 
constant g"^ — tt/IOO and mass scale ji/m — Q.l. 
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